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In this paper we analyze the interaction of a uniformly accelerated detector with a quantum 
field in (3+l)D spacetime, aiming at the issue of how kinematics can render vacuum fluctuations 
the appearance of thermal radiance in the detector (Unruh effect) and how they engender flux of 
radiation for observers afar. Two basic questions are addressed in this study: a) How are vacuum 
fluctuations related to the emitted radiation? b) Is there emitted radiation with energy flux in 
the Unruh effect? We adopt a method which places the detector and the field on an equal footing 
and derive the two-point correlation functions of the detector and of the field separately with full 
account of their interplay. From the exact solutions, we are able to study the complete process from 
the initial transient to the final steady state, keeping track of all activities they engage in and the 
physical effects manifested. We derive a quantum radiation formula for a Minkowski observer. We 
find that there does exist a positive radiated power of quantum nature emitted by the detector, 
with a hint of certain features of the Unruh effect. We further verify that the total energy of the 
dressed detector and a part of the radiated energy from the detector is conserved. However, this 
part of the radiation ceases in steady state. So the hint of the Unruh effect in radiated power 
is actually not directly from the energy flux that the detector experiences in Unruh effect. Since 
all the relevant quantum and statistical information about the detector (atom) and the field can 
be obtained from the results presented here, they are expected to be useful, when appropriately 
generalized, for addressing issues of quantum information processing in atomic and optical systems, 
such as quantum decoherence, entanglement and teleportation. 

PACS numbers: 04.62.+V, 04.70.Dy, 12.20.-m 



I. INTRODUCTION 



Inasmuch as studies of the interaction between a particle and a quantum field are basic to particle physics and field 
theory, understanding the interaction between an atom and a quantum field is essential to current atomic and optical 
physics research Q, Q, IE ^ IS- The interaction of an accelerated charge or detector (an object with some internal 
degrees of freedom such as an atom or harmonic oscillator) in a quantum field is a simple yet fundamental problem 
with many implications in quantum field theory |^ Q , thermodynamics |3 and applications in radiation theory 
and atomic-optical physics. 

It is common knowledge that accelerating charges give rise to radiation [lo| . But it is not entirely straightforward 
to derive the radiation formula from quantum field theory. How are vacuum fluctuations related to the emitted 
radiation? When an atom or detector moves at constant acceleration, according to Unruh ^2|> it would experience a 
thermal bath at temperature Tjj ~ fia/ {2'iTckB)^ where a is the proper acceleration. Is there emitted radiation with 
an energy flux in the Unruh effect? 

Unruh effect, and the related effect for moving mirrors studied by Davies and Fulling |0, were intended originally to 
mimic Hawking radiation from black holes. Because of this connection, for some time now there has been a speculation 
that there is real radiation emitted from a uniformly accelerated detector (UAD) under steady state conditions (i.e., 
for atoms which have been uniformly accelerated for a time sufficiently long that transient effects have died out), 
not unlike that of an accelerating charge 0, ■ In hght of pending experiments both for electrons in accelerators 
[T^ ITsl . IT^ and for accelerated atoms in optical cavities ^3 this speculation has acquired some realistic signiflcance. 
There is a need for more detailed analysis for both the uniform acceleration of charges or detectors and for transient 
motions because the latter can produce radiation and as explained below, sudden changes in the dynamics can also 
produce emitted radiation with thermal characteristics. 
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After Unruh and Wald's 'iq earlier explication of what a Minkowski observer sees, Grove [l9j questioned whether 
an accelerated atom actually emits radiated energy. Raine, Sciama and Grove pol | (RSG) analyzed what an inertial 
observer placed in the forward light cone of the accelerating detector would measure and concluded that the oscillator 
does not radiate. Unruh in an independent calculation, basically concurred with the findings of RSG but he also 
showed the existence of extra terms in the two-point function of the field which would contribute to the excitation of 
a detector placed in the forward light cone. Massar, Parantani and Brout [l^l (MPB) pointed out that the missing 
terms in RSG contribute to a "polarization cloud" around the accelerating detector. For a review of earlier work 
on accelerated detectors, see e.g., ^3J. For work after that, see, e.g., Hinterleitner [23|, Audretsch, Miiller and 
Holzmann [23 , Massar and Parantani |2g . Our present work follows the vein of Raval, Hu, Anglin (RHA) and Koks 
[27L I28L l2^ l30l | on the minimal coupling model and uses some results of Lin on the Unruh-DeWitt model [l^ls^ . 

With regard to the question "Is there a radiation flux emitted from an Unruh detector?" the findings of RSG, Unruh, 
MPB, RHA and others show that, at least in (1+1) dimension model calculations, there is no emitted radiation from a 
linear uniformly accelerated oscillator under equilibrium conditions, even though, as found before, that there exists a 
polarization cloud around it. Hu and Johnson emphasized the difference between an equilibrium condition (steady 
state uniform acceleration) where there is no emitted radiation, and nonequilibrium conditions where there could be 
radiation emitted. Nonequilibrium conditions arise for non-uniformly accelerated atoms (for an example of finite 
time acceleration, see Raval, Hu and Koks (RHK) |30||'). or during the initial transient time for an atom approaching 
uniform acceleration, when its internal states have not yet reached equilibrium through interaction with the field. Hu 
and Raval (HR) 27, 28] presented a more complete analysis of the two-point function, calculated for two points lying 
in arbitrary regions of Minkowski space. This generalizes the results of MPB in that there is no restriction for the 
two points to lie to the left of the accelerated oscillator trajectory. They show where the extra terms in the two-point 
function are which were ignored in the RSG analysis. More important to answering the theme question, they show 
that at least in (1-1-1) dimension the stress-energy tensor vanishes everywhere except on the horizon. This means that 
there is no net flux of radiation emitted from the uniformly accelerated oscillator in steady state in (H-l)D case. 

Most prior theoretical work on this topic was done in (1-1-1) dimensional spacetimes. However since most exper- 
imental proposals on the detection of Unruh effect are designed for the physical four dimensional spacetime, it is 
necessary to do a detailed analysis for (3+1) dimensions. Although tempting, one cannot assume that all (3+1) 
results are equivalent to those from (1+1) calculations. First, there are new divergences in the (3+1) case to deal 
with. Second, the structure of the retarded fleld in (3+1) dimensional spacetime is much richer: it consists of a bound 
field (generalized Coulomb field) and a radiation field with a variety of multipole structure, while the (1+1) case has 
only the radiation field in a different form. Third, an earlier work of one of us [sj] showed that there is some constant 
negative monopole radiation emitted from a detector initially in the ground state and uniformly accelerated in (3+l)D 
Minkowski space, and claimed that this signal could be an evidence of the Unruh effect. This contradicts the results 
reported by HR [23 and others from the (1+1)D calculations. We need to clarify this discrepancy and determine the 
cause of it, by studying the complete process from transient to steady state. In particular, since radiation only exists 
under nonequilibrium conditions in the (1+1) case, it is crucial to understand the transient effects in the (3+1) case 
to gauge our expectation of what could be, against what would be, observed in laboratories. 

In conceptual terms, one is tempted to invoke stationarity and thermality conditions for the description of an UAD. 
This is indeed a simple and powerful way to understand its physics if the detector undergoes uniform acceleration and 
interacts with the field all throughout (e.g., because of the stationarity of the problem in the Rindler proper time it 
is guaranteed that the total boost energy-operator is conserved). However, this argument is inapplicable for transient 
epochs during which the physics is quite different (see, e.g., the inertial to uniform acceleration motion treated in [3(|). 
Likewise, one can invoke the thermality condition (i.e., the thermal radiance experienced by a UAD is equivalent to 
that of an inertial detector in a thermal bath) to obtain results based on simple reasonings. But then we note that 
the thermality condition does not uniquely arise from uniform acceleration conditions. For example if the motion is 
rapidly altered j34i] the radiation produced can be approximately thermal. This thermality in emitted radiation (e.g., 
from sudden injection of atoms into a cavity) is similar to those encountered in cosmological particle creation |35| . 
but has a different physical origin from Unruh effect which is similar to particle creation from black holes (Hawking 
effect) (see, e.g., llllSa]). 

In terms of methodology, instead of u sing the more sophisticated influence functional method as in the earlier series 
of papers on accelerated detectors |2^ l3Clj | and moving charges [s^, IS US i our work here follows more closely the 
work of HR who used the Heisenberg operator method to calculate the two-point function and the stress-energy tensor 
of a massless quantum scalar field. In our analysis based on the (3+l)D Unruh-DeWitt detector theory we found the 
full and exact dynamics of the detector and the field in terms of their Heisenberg operator evolution, thus making 
available the complete quantum and statistical information for this detector-field system, enabling us to address the 
interplay of thermal radiance in the detector, vacuum polarization cloud around the detector, quantum fluctuations 
and radiation, and emitted flux of classical radiation. 

The paper is organized as follows. In Sec. II wc introduce the Unruh-DeWitt detector theory. Then in Sec. Ill 
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we describe the quantum dynamics of the detector-field system in the Heisenberg picture, yielding the expectation 
values of the detector two-point function with respect to the Minkowski vacuum and a detector coherent state in 
Sec. IV. With these results we derive the two-point function of the quantum field and describe what constitutes the 
"vacuum polarization" around the detector in Sec.V. Then in Sec. VI we calculate the quantum expectation values of 
the stress-energy tensor induced by the uniformly accelerated detector. This allows us to explore the conservation 
law and derive the quantum radiation formula. A comparison with the results in Ref. 31] follows in Sec. VII. Finally, 
we summarize our findings in Sec. VIII. 

II. UNRUH-DEWITT DETECTOR THEORY 

The total action of the detector-field system is given by 

S = SQ + S.i + 5/, (1) 

where Q is the internal degree of freedom of the detector, assumed to be a harmonic oscillator with mass mg and a 
(bare) natural frequency fio- 



(drQf-nlQ^ . (2) 



Here r is the detector's proper time. Henceforth we will use an overdot on Q to denote dQ{T)/dT. The scalar field <I> 
is assumed to be massless, 

S^ = - I d^x^d^'^'d''^. (3) 
The interaction action Sj for the Unruh-DcWitt (UD) detector theory has the form [l^ Is^. 

Si^Xo [ dr [ d^xQ{T)<i>{x)S'^ {x" - z'^(r)) . (4) 



where Aq is the coupling constant. This can be regarded as a simplified version of an atom. 

Below we will consider the UD detector moving in a prescribed trajectory z'^{t) in a four-dimensional Minkowski 
spacetime with metric 77^,^ = diag(— 1, 1, 1, 1) and line element ds^ = rip^udx^dx'^ . By "prescribed" we mean the 
trajectory of the detector is not considered as a dynamical variable, thus we ignore the backreaction effect of the field 
on the trajectory. (See Ref-js^] for an example where the trajectory and the field are determined self-consistently by 
each other.) The detector is made (by the act of an external agent) to go along the worldline 

z^(t) = (a""'^ sinhar, a""'^ coshar, 0, 0), (5) 

parametrized by its proper time r. This is the trajectory of a uniformly accelerated detector situated in Rindler wedge 
R (the portion i — a;^ < and t + x^ > of Minkowski space; see Chapter 4 of Ref.f^). 

III. QUANTUM THEORY IN HEISENBERG PICTURE 

The conjugate momenta (P(t), H(a;)) of dynamical variables (Q(T),$(a;)) are defined by 

5Q{t) 

By treating the above dynamical variables as operators and introducing the equal time commutation relations, 

[Q{t),P{t)] = in, (8) 

[$(t,x),n(t,x')] = ihd^{x-x'), (9) 
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one can write down the Heisenberg equations of motion for the operators and obtains 



d^Mr)+^lQ{r) = ^$(r,z(r)), 



{d^ - V2) <i>(x) = Ao / dTQ{T)S\x - z(t)), 



(10) 
(11) 

which have the same form as the classical Euler-Lagrange equations. 

Suppose the system is prepared before r = 0, and the coupling Si is turned on precisely at the moment r = 
when we allow all the dynamical variables to begin to interact and evolve under the influence of each other. (The 
consequences of this sudden switch-on and the assumption of a factorizable initial state for the combined system a 
quantum Brownian oscillator plus oscillator bath is described in some details in |4(l||\ By virtue of the linear coupling 
the time evolution of $(x) is simply a linear transformation in the phase space spanned by the orthonormal basis 
($(x), n(x), Q, P), that is, $(a;) can be expressed in the form 

Mt, x) = 1 d^x' [/*(t, X, x')l'(0, x') + f^it, X, x')n(0, x')] + /'^(x)Q(O) + /^(.t)F(O). (12) 

Here /*(a;,x'), f^{x,x'), f'^{x) and f^{x) are c-number functions of spacetime. Similarly, the operator Q{t) can be 
written as 



d^x' 



g*(r, x')$(0, x') + q^ir, x')n(0, x') + <z«(r)Q(0) + <z^(r)P(0), 



(13) 



with c-number functions q'^{T),q^{T),q'^{T,x') and g^(T, x'). 

For the case with initial operaters being the free field operators, namely, $(0,x) ~ $o(x), <I'(0,x) = no(x), 
Q(0) = Qo and P(0) — Pq, one can go further by introducing the complex operators 6k and a: 



$o(x) 

no(x) 



d^k 
(27r)3 V 2lj 

d^k fn 



'6k + e 



— ik-xf,t 



(27r)3 V 2lj 



k 



with 



Ikl, and 



(14) 
(15) 

(16) 



Note that, instead of fio, we use the renormalizd natural frequency Qr (to be defined in (|40|l ') in the definition of d. 
Then the commutation relations © and ^ give 



and the expressions H12|l and H13|l can be re-written as 

l'(t,x) = ^b{x)+^a{x), 
Q{t) = Qb{T)+QaiT) 



where 



Qbir) 

Qa{T) 



d^k 
(2^ 



251r™o 
d^k I h [ 



(27r)3 V 2uj 



/W(i,x;k)6k + /(-Ht,x;k)6t 
[r(i,x)a + r*(i,x)at] , 

gW(T,k)6k + <Z^"'(T,k)6l, 



2r2rTOO 



[q'^{T)d + q^*{T)d^ 



(17) 

(18) 
(19) 

(20) 
(21) 
(22) 
(23) 
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The whole problem therefore can be transformed to solving c- number functfons /''(a;) and (?''(t) from H10() and Hll() 
with suitable initial conditions. Since Q and $ are hermitian, one has /^"■' = (/'•"'''')* and g'") = Hence it is 

sufficient to solve the c-number functions /(+^(t,x;k), (7'^+)(r,k), /"(i,x) and ^"(t). To place this in a more general 
setting, let us perform a Lorentz transformation shifting r = to t = tq, and define 

7/ = r - To. (24) 

This does not add any complication to our calculation. Now the coupling between the detector and the field would 
be turned on at t — tq. We are looking for solutions with the initial conditions 

/(+)(i(ro),x;k) -e*-, (t(To), x; k) - -zc.e"^-, (tq; k) = (to; k) = 0, (25) 

r(t(To),x)=9ir(i(ro),x)-0, g''(To) = l, r{ro) ^ ~i^r- (26) 



A. Solving for and g*"*"^ 

The method to solve for / and q are analogous to what we did in classical field theory. First, we find an expression 
relating the harmonic oscillator and the field amplitude right at the detector. Then substituting this relation into the 
equation of motion for the oscillator, we obtain a complete equation of motion for q with full information of the field. 
Last, we solve this complete equation of motion for q, and from its solution determine the field / consistently. 

Ea. l(TT|) implies that 



(52-V2)/(+)(x;k)=Ao / dTS\x~z{T))q^+\r;k). (27) 
The general solution for /*^+) reads 



where 



is the free field solution, and 



f^+\x;k) ^ f^+\x;k) + fi+\x;k), (28) 
f^+\x;k) = 6^"^'+'^-'' (29) 

poo 

/+\x;k) = Xo dTGrct(x;z(T))g(+)(T;k) (30) 

J To 

is the retarded solution, which looks like the retarded field in classical field theory. Here w = |k| and the retarded 
Green's function Gret m Minkowski space is given by 

G,,t{x,x') = ^S{a)e{t-t') (31) 

with a = ~{x^ — x'^){x'^ — x'^)/2. Applying the exphcit form of the retarded Green's function, one can go further to 
write 

/i"^(x;k) = M^,W(,_;k), (32) 

where 



X = y/i-UV + p2 + a-2)2 + 4a-mV, (33) 

r- ^ -^\n^^{X-UV + p' + a-^), (34) 

V- = T- - To, (35) 

with p = \/ + 2^3 U = t — x^ and V = t + x^ . 

The formal retarded solution (|32|) is singular on the trajectory of the detector. To deal with the singularity, note 
that the UD detector here is a quantum mechanical object, and the detector number would always be one. This 
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means that at the energy threshold of detector creations, there is a natural cutoff on frequency, which sets an upper 
bound on the resolution to be explored in our theory. Thus it is justified to assume here that the detector has a finite 
extent 0(A^^), which will introduce the back reaction on the detector. 

Let us regularize the retarded Green's function by invoking the essence of effective field theory: 

GU^, x') = l.^K^e-^^'<^^e{t - t'). (36) 



(For more details on this regularization scheme, see Refs.|37l l39l|.) With this, right on the trajectory, the retarded 
solution for large A is 

/l+)(z(T);k) = ^ [ACgW(r;k)-a.g(+)(r;k)+0(A-i)] , (37) 

where C, = 2^/''r(5/4)/-y/7r. Substituting the above expansion into H10|l and neglecting 0(A^^) terms, one obtains the 
equation of motion for q^^^ with back reaction, 

(a2 + 279. + f^2)^(+)(^;k) = ^/(+)(z(r);k). (38) 

Too 

Fortunately, there is no higher derivatives of q present in the above equation of motion. Now q^~^^ behaves like a 
damped harmonic oscillator driven by the vacuum fluctuations of the scalar field, with the damping constant 



(39) 



Stttoq 

and the renormalized natural frequency 

nl^nl^f^. (40) 

In H38|l . the solution for (7'+^ compatible with the initial conditions g'^+-'(ro;k) = g(+^(To;k) = is 

g(+)(r;k) = ^ >] / dr'c,e-^(^-^')/i+^(z(r');k), (41) 



^ /'dr'c,e-^(^-^')/r'(^(r');k), 
where /q^^ has been given in (|29fl . c± and w± are defined as 



c±=±:-^, w±=-7±z!^, (42) 



with 



n = ^ni - 72. (43) 

Throughout this paper we consider only the under-damped case with 7^ < f2^, so f2 is always real. 

B. Solving for and 

Similarly, from H1(J|I . Hll|l. H18|l and (|19|l . the equations of motion for /" and q"^ read 

(92 _ v2) r(x) = Ao 1 dT5\x - z{T))q^{T), (44) 

(d'^ + nl)q^ir)^^nz{T)). (45) 
Too 

The general solutions for similar to (|28|l . is 

poo 

r(x) =/o"(x) + Ao / drG,ct(x;z(T))g''(r_) (46) 

To 
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However, according to the initial condition (|26|l . one has /q = 0, hence 



(47) 



Again, the value of is singular right at the position of the detector. Performing the same regularization as those 
for q^+\ (113) becomes (cf. (jSHJ) 



(48) 



which describes a damped harmonic oscillator free of driving force. The solution consistent with the initial condition 
q'^(To) = 1 and q''{To) = —iilr reads 



1 - 



n 



1 



Or + i"f 



^ — iQri 



(49) 



IV. TWO-POINT FUNCTIONS OF THE DETECTOR 



As shown in the previous section, as Q evolves, some non-zero terms proportional to $ and II will be generated. 
Suppose the detector is initially prepared in a state that can be factorized into the quantum state | q ) for Q and the 
Minkowski vacuum | Oj\/) for the scalar field that is, 



\ro)^\q)\ Om) 
then the two-point function of Q will split into two parts, 



(50) 



Q{r)Q{r') > 



OmI ( q I [Qb{T) + QaiT)\ [Qbir) + Q,(t)J | q ) |0m 
q\q){ Q{T)Q{r') + ( Q{t)Q{t') )^ (Om|Om) • 



where, from H19|l . 



Q{r)Q{T') = (Om \Qb{r)Qb{r')\ Om) 
Q{t)Q{t')), = {q\Qa{r)Qa{T)\q). 



(51) 



(52) 
(53) 



Similar splitting happens for every two-point function of $(a;) as well as for the stress-energy tensor. 

Observe that ( Q{t)Q{t') )^ depends on the initial state of the field, or the Minkowski vacuum, while ( Q{t)Q{t') )^ 
depends on the initial state of the detector only. One can thus interpret ( Q{t)Q{t') )^ as accounting for the response 
to the vacuum fluctuations, while ( Q{t)Q(t') )^ corresponds to the intrinsic quantum fluctuations in the detector. 

In the following, we will demonstrate the explicit forms of some two-point functions we have obtained and analyze 
their behavior. To distinguish the quantum or classical natures of these quantities, the initial quantum state | q ) will 
be taken to be the coherent state 



q) =e 



E 

Tl = 



(54) 



where | n ) is the n-th excited state for the free detector, and a = qQyJVlr/2h with a constant go- The representation 
of I g ) in Q-space reads 



O, 



1/4 



(55) 



which is a wave-packet centered at go with the spread identical to the one for the ground state. 



8 



A. Expectation value of the detector two-point function with respect to the Minkowski vacuum 

Along the trajectory z^{t) in Q, performing a Fourier transformation with respect to r on ()29|l. one has 



where the frequency spectrum of the Minkowski mode from the viewpoint of the UAD, 



V5(K,k) 



-TTK/2a 



LO — ki 



K 



na \uj + ki 

is not trivial any more. Given the result of the integration, 

f h(fk / , X / , X ^ 
J (2^^("'^)^ = (2^2i_,-2W.- 

a Planck factor with the Unruh temperature a/27r emerges. Then from H52f) . (|41|l and (|56|l . one has 

d^k 



QMQiv') 



(27r)32tj 



gW(T;kV-)(r;k) 



(-)( 



(2 



nyml^J 1 



(56) 



(57) 



(58) 



(59) 



(60) 



where the integrand has poles at k = ±fi — 17 and k — ±ina, n £ N . Let Tq < tq < t < r' and taking the coincidence 
limit, one obtains 



Qiv? 



lim l({Q(r;),Q(V)}), 

\2 

^0(,7)Rc {(Ao - lna)e-27. gj^^ 



^^7+^0(6-'''') / in 



1 

in 

7 



7 + if2 + a \ 7 



if] 



F-7-»n(e''"') 
7 + zfi — a 



1 -e 



-2jf2?7 



iTT coth — (ri — 17) 
a 



Here Fs (t) is defined by the hyper-geometric function as 



Fs{t) = 2Fi(l + -,l,2+-;t] , 
a a 



(61) 



(62) 



and 



V". = ^ (1 + ^) 



(63) 

is the poly-gamma function. The divergent Ap-term is produced by the coincidence limit: as 77' — > 77, Aq ^ — 7e — 
In Itq — To I with the Euler's constant 7e. Since Itq — ro| characterizes the time scale that the interaction is turned on, 
Aq could be finite in real processes. In any case, for every finite value of Aq, the first line of the result in 161|l vanishes 
as 77 ^ 00. 

In FIG.n we show the ( Q{ri)^ )^ without Ao-term in dotted line. Roughly speaking the curve saturates exponen- 
tially in the detector's proper time. As — > 00, { Q{ri)^ )y saturates to the value 



lim ( Qiijf )^ = -Re 



-f + ifl 



(64) 
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FIG. 1: Plot of ( Qinf >v (dotted line, Eq.||nU with Ao-term excluded), ( Q{r]f )^'" (dashed line, Eq.^), and the sum of 
these two (( AQ{r])^ ), solid line). Here we have taken a = 1, 7 = 0.1, fl — 2.3 and mo = 1. 



For 7 < a, the time scale of the rise is about I/27, which can be read off from the e '^^^ in (|61|l . From there one can 
also see that the small oscillation around the rising curve has a frequency of 0{fl). 

For ( Q{'rj)Q{rj) )^, it will be clear that what is interesting for the calculation of the flux is the combined quantities 
like ( Q{ri)Q{ri) )^ + ( Qiv)Qiv) )v Notice that 



(65) 



With the result of ( Q{riY this calculation is straightforward. Let us turn to the two-point functions of Q. Similar 
to ( Q{ri)^ )y, one has 



(27r)32cj 

E 



(T;k)g(-)(T;k) 
kcIk 



1 



27rK/a (^^^ _|_ in){w*, - in) 



W 



* w',(t' —t!,) 

Wj,e °' 



me 



iK(T'-Ta) 



(66) 



(67) 



from (|41|) . H56|) and H58|l . The coincidence limit of the above two-point function reads 



(27rmori)' 



r6l(77)Re | (Ai - In a)^'^ + (Ao - In a)e"2'^'' [VL cos 1777-7 sin VL-qf 



+ ^(7 + zr!)2e-(^+'^)'' 



m 

1 



7+iO 



(in 






he" 


\~' 




-277, 


fiVt 




\ 7 



^ + iO, + a 
ifl 



7 — ^0 



7 



7 -1- — a 



^ in 



JflTJ 



^-2.n^^^ z7rcoth^(f7-i7) | 



(68) 



where Ai = — 7e — lim^r'^r In I''" ~ ''"'I can be subtracted safely. This will be justified later. 

The subtracted ( Q{r])^ )^ is illustrated in FIG.[21(in which Ao-term has also been excluded). One can immediately 
recognize that ( Qii])^ )^ ~ Ini] when rj approaches zero; a new divergence occurs at = 0. Mathematically, this 
logarithmic divergence comes about because the divergences in the hyper-geometric functions in (|68|l do not cancel 
each other, unlike in ( Q{riY Physically, this divergence at the initial time tq could be another consequence of the 
sudden switch-on at r = tq or 77 = 0. We expect that these ill-behaviors at the start could be tamed if we turn on 
the coupling adiabatically. (See |40| for a discussion on this issue.) 

For large -q, the behavior of the dotted curve in FIG. |21 is quite similar to the one in FIG. Q] for ( Q{tY )v- It 
saturates to 



hm ( Q{T^f 
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27rmofi 
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FIG. 2: Plots of ( Q{-nf >^ (dotted line, Eq.JSHJ), ( Qivf )T (dashed line, Eq.JZHJ) and their sum ({ AQ{T^f ), solid line). 



Comparing H68|l and (|61|1 . their time scales of saturation (I/27 for 7 < a) and the frequency of the small ripples on 
the rising curve (0(il)) are also the same. Note that when 7^1 and a is finite, Ht)9|) and (|54|l implies 



(70) 



which justifies the subtraction of Ai-term in (|68|l . 



B. Expectation values of the detector two-point functions with respect to a coherent state 



We now derive the expectation values of the detector two-point functions with respect to the coherent state H54|l. 
Subsituting into itS^ and using and one finds that 



where 



( Q(r)Q(r') = ( Q(r)g(r') )r + ( Q{T)Q{r') , (71) 



[Q{r)Q{T'))r ^ :^^q''ir)q''*(r'), (72) 

( Q{r)Q{r') )t ^ ^Re [g'^(r)] Re [g«(r')] = Q{r)Q{T'), (73) 
Too 



with the mean value 



Q{t) = ( Q(t) ) = ^^^(,7)e-T" (cos r^T? + H sin flT]) . (74) 

^TOq V ll ' 

While the "qm" term is of purely quantum nature, the "cl" term is of classical nature: Q is real and ( Q{t)Q{t') 
does not involve h. Thus the "cl" term is identified as the semiclassical part of the two-point functions. The coincidence 
limits of the above two-point functions are 

( Qivf >r = J^^-''' [nl-j'cos2nr^ + jnsin2nrj] , (75) 



Similarly, it is easy to find ( Q{t)Q{t') )^ = ( Q{t)Q{t') )^'° + ( Q(t)Q(t') );;': 

{Q{r)Q{r'))T = 7^i\r)q'*{r'), (76) 

( Q{r)Q{r') = Q{r)Q{r'), (77) 

and their coincidence limit, 

( Qivr )r = [^r - 1' cos 217^ - ^nsm2nr^] , (78) 
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FIG. 3: The semiclassical part of the two-point funciton, { )|^' = (see Eg. 1)73^ and below). Its behavior is quite different 
from the quantum part shown in the previous figures. Here we take qo — 1, with other parameters unchanged. 
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FIG. 4: ( AQ(oo)2 ) = ( Q{oof >^ against the proper acceleration a (solid line, Eg. 1641 ') with other parameters the same as 
FIG. □ For small a, the value of ( Q{oof )^ is less than ( AQ(0)^ ) = ( Q{Of )^ (dotted line, see Eq.JTSJ)). For large a, 
{ Q{oo)^ )^ is nearly proportional to a. 



and ( Qirj) 



2 \cl 



Q(7?)2. Also one has ( Q{t)Q{t) )^ + ( Q{t)Q{t) = 9, ( Q{i 



Note that the above two-point functions with respect to the coherent state are independent of the proper acceleration 
a. ( Qit])^ )g[™ and the variance (squared uncertainty) of Q, 



[Qiv) - Qiv)? 



Qiv) 



2 \qm 



(79) 



have been shown in FIG.n ( Q{ri)^ decays exponentially due to the dissipation of the zero-point energy to the 
field. As ( Q{r])^ )g[™ decays, ( Q{ri)^ )^ grows and compensates the decrease, then saturates asymptotically. Similar 
behavior can be found in FIG. |2 in which 



2 \qm 



(80) 



is illustrated. In FIG.|31we show the semiclassical two-point funciton {Q'^Ya- Its behavior is quite different from the 
quantum part shown in the previous figures. 



C. Late-time variances and the proper acceleration 



The saturated value ( Q(oo)^ )^ in Ea. (|54|l is the late-time variance of Q, namely, ( AQ(oo)^ ) = ( (5(cx3)^ )^. Its 
dependence on the proper acceleration a is shown in FIG. 01 

One can see that, when a is large, ( Q(oo)^ )^ is nearly proportional to a, while in the zero-acceleration limit a — )■ 
with Tj ^ (27)"^ In I lna|, the saturated value goes to a positive number. From (|64|l and (|75|) . one finds that 



g(oo) 



lim , „ , , 

a-o ( g(o) 



2 



In 



(81) 
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thus ( AQ(cx))^ ) = ( (5(oo)^ )^ is smaller than ( AQ(0)^ ) = ( Q(0)^ for every 7 > when a —>■ 0. In other words, 
for a non-accelerated detector, whose Unruh temperature is zero, the variance of Q in the detector-field coupled 
system is still finite and smaller than the one for the ground state in the free theory. 

Actually, ( AQ(oo)^ ) will become smaller than ( AQ{0)^ ) whenever a is small enough. Observing FIG. 01 there 
is a critical value of a that gives the late-time variance identical to the initial one (a = ~ 3.0447 in FIG. 0J. 
Does this mean that the Q-componcnt of the final wave-packet with acr is in the original ground state of the free 
theory? The answer is no. What happens is that the quantum state of Q has been highly entangled with the quantum 
state of $ at late times, and the value of 1/ ( AQ{oo)'^ ) simply represents the width of the projection of the whole 
wave-packet (in the Q-^ representation of the state) onto the Q-axis. There is actually no factorizable Q-component 
of the wave-packet, and the final configuration of the wave-packet in Q-^ space looks totally different from the initial 
one. Indeed, with the same critical value of a, ( A(3(oo)^ ) — ( Q{oo)'^ )^ is not equal to ( AQ{0)^ ) = ( Q{0)^ )2™ for 
every 7^0. 

But one can still imagine that, at = 0, the coherent state for the free detector is an ensemble of particles with a 
distribution function like | ( Q | g ) P in Q-space. \/ ( QivY )a'" width of this distribution function. When 77 > 0, 

due to the dissipation which comes with the coupling, all particles in the ensemble are going to fall into the bottom 
of the potential of Q, so ( Qirff' )^"^ shrinks to zero. On the other hand, the vacuum fluctuations of the field act like 
a pressure which can push the ensemble of particles outwards, so that the width of the projection of the wave-packet 
in Q-space, { AQ^ ), remains finite. A larger a gives a higher Unruh temperature, and a higher outward pressure, 
so eventually the wave-packet reaches equilibrium with a wider projection in the potential well of Q. 



D. Shift of the ground state energy 



A natural definition of the energy of the dressed detector (a similar concept is that of a "dressed atom", see e.g., 
Ref. ^,^) is 



mo 



(82) 



with ( Q2(^) ) = 
can see that Q, 



{ QHv) )v + ( Q\v) )a and ( Q\r]) ) = ( Q\r^) + ( Q^{r]) )^ according to In FIGs. HE one 

Q^iv) )a'" and ( Q^{r}) )|^™ eventually die out. So the late-time energy of the dressed detector is 



E{oo) 



mo 
2 
fi 

2^ 



{a - 2Re [(7 + - 27 In a} 



(83) 



from and This is actually the true ground-state energy of the dressed detector, with the vacuum fluctuations 
of the fleld incorporated. The first term in E{oo) could be interpreted as the total energy of a harmonic oscillator in 
thermal bath, kBTjj, with the Unruh temperature Tu — Tia/2'KkB- 

The ground-state energy of the dressed detector is not identical to the one for the free detector, Eq — 'hQ.r/2. In 
particular, if a is small enough, the subtracted E{oo) is lower than £^0, though there is an ambiguity of a constant in 
determining the value of the energy. This is analogous to the Lamb shift in atomic physics 0, |22| . 



V. TWO-POINT FUNCTIONS OF THE QUANTUM FIELD 

Similar to the two-point functions of the detector, for the initial quantum state ()50fl , the two-point function of $ 
could be split into two parts, 

( $(X)$(X') ) = (Om \{q\ [^a{x) + $6(X)] [$a(x') + ^x' )] \q)\ Om) 

= G^{x,x') + G^{x,x'), (84) 

where, from l(T^ and (fTH|l . 

^^^f(+\x;k)f(-\x';k), (85) 

Ga(a;,a;') ^ (Om | Om) ( g | <i>6(:r)$,(a;') U ) = ^-^/^ (x)r*(x'). (86) 
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Eas. H28|l - H30|l and H47|) suggest that Gv accounts for the back reaction of the vacuum fluctuations of the scalar field 
on the field itself, while Ga corresponds to the dissipation of the zero-point energy of the internal degree of freedom 
of the detector. 

Substituting ^ into (|H3J), Gv can be decomposed into four pieces, 

G^{x, x') = G™(a;, x') + G^i(a;, x') + Gl°{x, x') + Gl\x, x'), (87) 

in which Gl^ are defined by 

h(fk 



Gv (x^x ) 



(27r)32w 



f^+\x;k)fl-\x',k), 



with i,j — 0,1. G^° is actually the Green's function for free fields, which should be subtracted to obtain the 
renormalized Green's function for the interacting theory, namely. 



G„„(x,a;') = ( Hx)^{x') ) -G^°(x,x'). 



(89) 



Since G'^^{x,x') = [G^'^ {x' , x)]* by definition, it is sufficient to calculate Gv^(a;,a;') and G]^^{x,x') in the following. 
The structure of G^^ is quite simple. Comparing H32I) . H59|l and the definition (|88|l . one concludes that 



Gv {x, x') 



X2 
-^0 



- ( Q{r-)Q{rL) 



(90) 



{2TT)^a^XX' 

The result can be substituted directly to get the coincidence limit of G^^ . 

By definition, G^° accounts for the interference between the retarded solution /l^-* and the free solution /o^^ Since 

we are interested in the coincidence limit of G^'^, and /j;^'' vanishes in the L- wedge {U > 0, V < 0) and P-wedge 
([/, < 0) of Minkowski space, below only the Gv''(a;, x') with x and x' in the F- wedge {U, V > 0) and R- wedge would 
be calculated. 

It has been given in Ref.|23 that 



(27r)32w 



(^(K,k)e' 



ot—i\z.-yi 



e*'"^dT/(27r) 



(xi - z^{T)y + p^-{t- zO(t) + le) 



2-KaX (1 - e-2^«/'») 



Z(K)e*'' 



(91) 



where e > 0, Z{k) — 1 and e '^'^Z" for x in R and F-wedges respectively, X and r_ were defined in Ip?!^ and 
and 



l.ln.« 



a 2\U\ 



Hence, from (I32|l and H29|) . one has 

M20(r7_) 



Gi"(x,a;') 



(27r)3TOoa2XX' 



(92) 



(93) 



with r/± (x) = r-t (x) — tq . The coincidence limit of G}^ reads 

Gi°(x,x) ^ hm i(Gi°(x,x') + Gi°(x',x)) 
?iA20(r;_ 



(27r)3mor2a2X2 



Re <^ itjj^+in + 



J + iQ + a 



e"(^+*"+'^)''-F^+,o(e~"''-) 



-(±)e-(''+*")''--"''+F^+,n(±e"'"'+) ± e-''(''+-''-)i^-,+io(±e-''(''+-''-)) } 



(94) 



with + and — for x in R and F-wedges, respectively. Near the event horizon ?7 ^ 0, 77+ diverges, and the last two 
terms in (|94|l vanish. 

As for Ga(x,x'), since /q = 0, one has G^^ = G^^ = 0, and only G^^ contributes to Ga- Inserting (g7|l into (jSHJ 
and comparing with H71|l. one finds that 

Ga(x,x') 



-(Q(r-)Q(/) 



i2iT)^a^XX"'°''''''°'''-' ^^^^ 
It can also be divided into a quantum part G1™{x, x') and a semiclassical part Ga (a;, x') according to H71|) and below. 
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FIG. 5: Plot of GJ," + Gv^ (dashed line), Gl^ + G^™ (dotted line) and their sum (solid line) against V near the event horizon 
i7 = 0. Other parameters are the same as those in FIG. One can see the feature that positive G" + Ga™ is screened by 
negative Gl° + G°\ 




FIG. 6: Plot of X^{Gl° + G°^) (dashed line), X^{Gi^ + GT) (dotted hue) and their sum (solid line) against U at V = 15 and 
p = 0. The cusp in the right (U « —0.07 in this plot) locates at the position of the detector with X^{GI'^ + G^^) — 0. It is due 
to the weaker divergence than for Gi° + G^^ . 



A. Effects due to the interfering term 

In our (3+1) dimensional UD detector theory, the coincidence limit of the quantum part of Gj-on reads 

^ G2'^{x,x) + Gl\x,x)+Gl°{x,x) + G°'^{x,x), (96) 

owing to (IH3I), 1113 and (jHU). Collecting the results in ^U^, and it is found that G^™(x,x) is singular at 
X y{T), and one has to be more cautious. 

As can be seen from I^HI and Gl^{x, x) and G2™{x, x) look like the squares of the retarded field with effective 
squared scalar charge ( Q{t)'^ )^ and ( Q{tY )a™j respectively. Since the detector is accelerating, these two terms 
do carry radiated energy (this will be shown explicitly later). The interfering term G]^{x,x) + G^^{x,x), is more 
intriguing: At first glance, it acts like a polarization in the medium, which screens the radiation field carried by 
G]^{x, x) and G^™(a;, x). However, the interfering term G^" + G°^ does not respond to G^™ at all - it is independent 
of /° and impervious to any information about the quantum state of Q. Hence the interfering term cannot be 
interpreted as the polarization in the medium. The total effect is simply a destructive interference between the field 
induced by the vacuum fluctuations, and the vacuum fluctuations themselves. For physical interpretations one should 
group G\^ + G°^ and G\^ together and leave G!J™ alone. 

These quantities, together with their sum, are illustrated in FIGs. |31and|ni In FIG.|31 one can see that, soon after 
the coupHng is turned on at y = 1/a, G^° + G°^ build up and pull the sohd curve down. Observing the time 
scale (in proper time of the detector) of this pull-down is about 1/(7 + a), which is shorter than the time scale I/27 
for ( Q{tY )v and ( Q{tY , since we take 7 = 0.1 < a — 1 here. In FIG. El one can also see that G^™ + G" 
diverge as around the trajectory, while the divergence of G}^ + G^^ as AT ^ is a bit weaker than A~^, such 
that A^(Gy° + G"^) goes to zero on the trajectory of the detector. 
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B. What exactly is the "vacuum polarization cloud" around the detector? 

In prior work for (1+1)D spacetime the counterpart of Gren(x,x) has been considered as evidence for the 
existence of a "vacuum polarization cloud" around the detector [23,|22,|23- This is because Gren{x,x) around the 
detector does not vanish even after the system reaches equilibrium, it exchanges particles with the detector, and 
the mean energy it carries is zero. Nevertheless, vacuum polarization is a concept pertinent to field-field quantum 
interacting systems. In quantum electrodynamics, electrons are described in terms of a field, which distributes in the 
whole spacetime, so vacuum polarization is pictured as the creation and annihilation of virtual electron-positron pairs 
everywhere in spacetime. These virtual electron-positron pairs do modify the field strength around the location of 
a point charge, yielding a non-vanishing variance of the electromagnetic (EM) field. But in the UD detector theory, 
at the level of precision explored here, the detector-field interaction (hence the virtual processes) only occurs on the 
trajectory of the detector. There is no virtual detector or scalar charge at any spatial point off the location of the UD 
detector. 

Hence in UD detector theory "vacuum polarization cloud" is not a precise description of G,:cn{x,x) in steady state. 
At late times G,-enix, x) simply shows the characteristics of the field in the true vacuum state, in contrast to ( Q(oo)^ ) 
for the true ground state of the detector. 



VI. RADIATED POWER 

In classical theory, the modified stress-energy tensor for a massless scalar field $ in Minkowski space is 



(97) 



where f is a field coupling parameter, set to zero here. Denote v^^ = dz'^/dr as the four velocity of the detector, and 
define the null distance r and the spacelike unit vector u'^ by — z^(r_) = r{u'^ + v^{t-)) |41| with normalization 
u^u'' = 1 and v^uf^ — (see FIG. [7||. Then the stress-energy tensor for the classical retarded field <l>rGt induced by 
the UD detector moving along the trajectory z''(r_) can be written as [sif 

1 



Tf^A^rctix)](=0 
1 



(4^) 



1 r 

2 



Q{t-) + Q{T-)apuP {Uf_, + Vf^){u^ +v^) 



(98) 



The 0{r ^) term in the above expression corresponds to the radiation field, which carries radiated power given by 



dW 



rad 



dr- 



lim / r^d^liiu''T^,^v"(T^) 



(47r)2 



dVL, 



Q{T-) + Q{T^)apuP 



(99) 



to the null infinity of Minkowski space. Here the Q-term corresponds to dipole radiation (? = 1, m = in multipole 
expansion of the radiation field) with the angular distribution apU^ = a cos 9 [44| . while the Q-term corresponds to 
monopole radiation {I — 0) isotropic in the rest frame instantaneously for the UD detector at r_. The solid angle 
dQu could be further integrated out, then one obtains the classical radiation formula 



dW 



rad 



x2 r 

47r 



(100) 



which is the counterpart of the Larmor formula for EM radiation. The second term is the usual radiation formula for 
the massless scalar field emitted by a constant, point-like scalar charge in acceleration [i^ . 

Naively, one may expect that the quantum version of the radiation formula could look like (Ao/47r)[( (5^(t_) ) -|- 
(a^/3) ( Q^{t-) )]. In the following, we shall calculate the quantum expectation value of the flux T^^, from which we 
will see that the quantum radiation formula is more complicated than expected. 
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FIG. 7: Definitions of p, 6 and r. 



A. expectation value of the stress-energy tensor 



The expectation value of the renormaUzed stress-energy tensor (T^^)^,^^ is obtained by calculating 



d d 1 d d 



Grcn ^ ) 1 



(101) 



according to with C = 0- With the results in the previous section, it is straightforward to obtain (T)„y)j.(,,j induced 
by the UAD: 



(V['i'(^)])rc„ = 



(27r)2a2x2 



X 



27rmo 



(102) 



Upon collecting IjAip and HA2p as well as those from C^. Here ( = ( ) + (?i/27rmo)9 , ( )j 



Q2 ) + (n/2^mo)exjf and ( QQ 



QQ ) + (?i/27rmo)9x- with 6^ defined in HASp - ljATp . To sec the properties 



of quantum nature, we define the total variances by subtracting the semiclassical part from ( • • • )^^^ 

hOxxiT-) 



as 



^(r_) = (AQ^(r_)> + 



( AQ'(r-) 



Q"(r_) = (AQ2(^_ 



27rmo 
27rmo 



(103) 
(104) 



Their evolution against ?]_ arc illustrated in FIGs. |Hland|H| 

In our case, the Minkowski coordinate (U, V, p) of a spacetime point in F and R-wedge can be transformed to the 
coordinate (r, t_ , 6) by 



p = r sin W, 

V = re"''- [1 + C0S6I+ (ar)-i] , 
U = re-"''- [1 -cos 61- (ar)-i] , 



(105) 
(106) 
(107) 
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FIG. 8: The total variance ( AQ^ )^^^ (Ea. ^Wd^ ) near the event horizon for the detector (77+ 00). This plot is virtually the 
same as FIG. [^except the "time" variable is rj- here. The values of parameters are still the same as before. The total variance 
finally saturates to the value ha/2-KmoQ.'^. One can compare this plot with FIG.Qdirectly and see the suppression. 



SO that X — 2r/a. Also one has 

u^{t^) = I r COS sinh aT_ , r cos 6* cosh aT_ , — , — j (108) 



with p — ^y(x^Y+{x^ . Now Ea. ()lfl2|l can be directly compared with (|98|l and (|99|l . One can see clearly that the 
( Q{t-)^ )^^^ term has the same angular distribution as the one for the term in H98() . hence would be recognized 
as a monopole radiation by the Minkowski observer. The angular distributions of the remaining terms in H1U2|) are, 
however, much more complicated because of their dependence on ?7+(r', r_, 0). 



B. Screening 



We have mentioned in the previous section that G^J™ and Gl} in (|95|l and H9U|) carry radiated energy, now this 
becomes clear. Observing that what correspond to 9^9,^/ [G^™(x, x') + Gl^{x,x')] are those proportional to ( • • • ) in 
( • • • terms of (|1U2|I . These terms contribute a positive flux. Nevertheless, due to the presence of the interfering 
terms Qij, most of this positive flux of quantum nature will be screened when the system reaches steady state as 
ry_ ^ 00. 

As shown in FIG.|H1 the total variance ( AQ^ )^^^ near the event horizon U — drops exponentially in proper time 
(power-law in the Minkowski time) after the coupling is turned on. Note that ( A(5^(?7_(x)) )^^^ is proportional to 
GJ?™(x,a;) defined in 19t)|) . and X is independent of V on the event horizon, so FIG.|Hlis virtually the same plot as 
FIG. iniexcept that the time variable here is rj-. Thus, similar to the behavior of G'}^^{x,x) near the event horizon, 
Qxx Gj" -|-G^^) builds up and the total variance ( A(5(r_)^ )^^^ is pulled down during the time scale 1/(7 + 0) (for 
7 < a) according to 1)61(1 . H75|l and l(A4p . Then ( AQ(t_)^ )^^^ turns into a tail (77- > 1 in FIG.|Hl which exponentially 
approaches the saturated value ha/2TTmoill with the time scale I/27. 

For ( A(5(t_)^ and ( AQ(t_)A(5(t_) )^^^, their behaviors are similar (see FIG. ISJ. In particular, 

( A(3(t_)^ = ( AQ(t_)^ ) + (n/27rmo)6 goes to zero at late times from (jSHl, (CHI and (|X3|) with 77?_ > 1, so 

the corresponding monopole radiation vanishes after the transient. 

From the calculations of Ref.jsJl based on perturbation theory, it was suggested that the existence of a monopole ra- 
diation could be an experimentally distinguishable evidence of the Unruh effect. Here we find from a non-perturbative 
calculation that, in fact, only the transient of it could be observed. (A comparison of both results will be given in 
Sec. VII.) This appears to agree with the claim that for a UAD in (l-t-l)D, emitted radiation is only associated with 
nonequilibriuni process |33|. The negative tail of ( AQ{t^ )^ )^^^^ in FIG. El and the corresponding quantum radiation 
could last for a long time with respect to the Minkowski observer i^^^ ^j^jg jg essentially a transient. 

The interference between the quantum radiation induced by the vacuum fluctuations and the vacuum fluctuations 
themselves totally screens the information about the Unruh effect in this part of the radiation. 
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FIG. 9: The total variance ( AQ^ >t^t (Eq.lESJ) with the same parameters. Note that ( AQ^ )^^^ is independent of Tj^, and 
this plot is not restricted around the event horizon. One can compare with FIG.|5|and see the suppression. 



C. Conservation between detector energy and radiation 



What is the physics behind the interfering term in ( Q^irj) 
can show that, 



Bin.) - Eirjf) = / drjiQ'ir,)) 



By inserting our results into H82|l and H1U4|) . one 



3. 


rvf 


Aw 




Ai 


rVf 


47r 





(109) 



for all proper time interval after the interaction is turned on (77/ > 77^ > 0). The left hand side of this equality is the 
energy-loss of the dressed detector from rji to 77/, while the right hand side is the radiated energy via the monopole 
radiation corresponding to ( Q^irj) during the same period. Therefore (|l()9|l is simply a statement of energy 

conservation between the detector and the field in this channel, and the interfering terms O must be included so 

that ( Q^{ri) )^^^ is present on the right hand side instead of the naively expected ( (5^(?7) ). A simpler but more 
general derivation of this relation is given in Appendix IbI Ea. (|109|) also justifies that (|82|l is indeed the correct form 
of the internal energy of the dressed detector. 

With the relation H109(l we can make two observations pertaining to results and procedures given before. First, 
while the Ao-terms in ( Q'^{t) )^ fEa. (|61|l l and ( Q^{t) )^ fEa. (|68|l l are not included in any figure of this paper, they 
are consistent with the conservation law H109I) . Actually the Ao-term in (|()1|) satisfy the driving- force-free equation of 
motion H48|l . just like the semiclassical Q does. 

Second, Ea. (|109p implies that all the internal energy of the dressed detector dissipates via a monopole radiation, 
and the external agent which drives the detector along the trajectory Q has no additional influence on this channel. 



D. Quantum radiation formula 

Transforming (|102|l to the form of (|98|l by applying (|105|) - (|107|l . one can calculate the radiation power 

'nuu>^{T^u),,,y{r-) (110) 
following a similar argument in classical theory. Before calculating, let us observe the behavior of the steady-state 
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T^i, )rcn the forward light cone. As r increses, the developments of two terms in late-time r^u^ ( T^^ 



are illustrated in FIGs. ^]and^2 It turns out that both are regular and non- vanishing at the null infinity of Minkowski 
space (r 00) even in steady state (777- 3> 1). FIGs. 1101 and ITTl also indicate that, near the the null infinity of 
Minkowski space, almost all the equal-r surface lies in the F-wedge, except the region around 6* = is still in the 
R- wedge. The contribution to the integral around 9 — can be totally neglected because the value oir'^u^ ( 



rcn 



V 



is regular there while the measure for this portion in the angular integral is zero when r — > cxd. So the radiation power 
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FIG. 10: Angular distributions of ( AQ'^ at late times (^Ea. p()8|l with 7r_ 3> 1) as the null distance r increases. The values 
of parameters are the same as before except that here we choose a = 2. The horizontal line indicates the saturated value of 
ha/2Tvmoflr shown in FIG.|H|. When r starts with 0, ( )^^^ decreases from ( Q^{oo) )^ (top-left). The light cone hits the 
event horizon at 6 — tv when r — l/{2a) = 1/4, and ( )^^^ always has the value at the event horizon (top-right, 

the curve and the horizontal line intersect right at the event horizon). As r further increases, { )^^^ sinks more and more 



(bottom- left), and some oscillations begin to develop near — n. Finally ( )t 
with the value smaller than ha/2-K'rrn)Q.^ whenever 9 ^ Q (bottom-right). 



is non- vanishing at the null infinity r oo, 



can be written as 
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by inserting H1U2|I . (|108|l and w^(t_) = (cosh ar_, sinh ar_, 0, 0) into (|110|l . This is the quantum radiation formula for 
the massless scalar field emitted by the UAD in (3-|-l)D spacetime. 

At late times, while ( )^ ^ ceases, it still remains a positive radiated power flow 
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(112) 



to the null infinity of Minkowski space. Thus we conclude that there exists a steady, positive radiated power of 
quantum nature emitted by the detector even when the detector is in steady state. 

For large a, the first term in (|112() dominates, and the radiated power is approximately 
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(113) 



where Tjj is the Unruh temperature. This could be interpreted as a hint of the Unruh effect. Note that it does not 
originate from the energy flux that the detector experiences in Unruh effect, since the internal energy of the dressed 
detector is conserved only in relation to the radiated energy of a monopole radiation corresponding to ( )tof 
Learning from the EM radiation emitted by a uniformly accelerated charge 0, ^| , we expect that the above non- 
vanishing radiated energy of quantum origin could be supplied by the external agent we introduced in the beginning to 
drive the motion of the detector. Further analysis on the quantum radiations of the detector involving the dynamics 
of the trajectory is still on-going. 
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FIG. 11: Angular distributions of x = u'^{'ri-V+ + V+V-)''^'^^+- ('Ea. l|A7|l ') at late times as r increases. Parameters are the 
same as those in FIG. IlUl When r starts with 0, x grows from (top- left). After the light cone hits the event horizon, x keeps 
regular on the event horizon (top-right, the vertical line indicates the position of the event horizon). As r further increases, a 
nodal point enters from the right (bottom- left). Then more and more nodal points can be seen. Finally x is non- vanishing at 
the null infinity r oo, with infinitely many nodal points close to ^ = tt (bottom-right). 



VII. COMPARISON WITH EARLIER RESULTS 



We can recover the results in Ref.jsJ, which is obtained by perturbation theory, as follows. In Ref.'si'l, the first 
order approximation of the flux ( T**^ )^^^ through the event horizon U = has been calculated. Here, the expectation 
value of r"^ near the event horizon reads, from (|102|l . 
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in which 0-| and Q+x terms vanish. By letting 7 — > with rj- finite, then taking rj^ 00, the total variance H104|l 

becomes 
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owing to lt)8|) . H78|l and (|A3|) . This is identical to the corresponding part of Eq.(66) in Ref . [sif . 
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(116) 



by noting that there, mo = 1, Y1,e I (-^0 \Q{^)\ E) P — {Eq |Q(0)^| Eq) = h/2£, and e there is equal to O^/a here. 

The monopole radiation corresponding to H116() looks like a constant negative flux since £ > 0. Accordingly it 
was concluded in Ref. j3ll| that such a quantum monopole radiation could be experimentally distinguishable from the 
bremsstrahlung of the detector. At first glance this constant monopole radiation seems to contradict the knowledge 
gained from (1-f 1)D results. But actually similar results for (l-l-l)D cases were also obtained by Massar and Parentani 
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(MP) [23, who found that a detector initially prepared in the ground state and coupled to a field under a smooth 
switching function does emit radiation during thermalization. They pointed out that the radiated flux in what they 
refered to as the "golden rule limit" (77 — > 00 with 'yn smalh while the switching function becomes nearly constant) is 
approximately a constant negative flux for all > . In spite of the long interaction time and the nearly constant 
radiated flux, the detector will remain in dis-equilibrium. 

Note that the initial conditions in are similar to those in Sec. II of MP, and the limiting condition for obtaining 
(|115|l is exactly what MP assumed there. Hence, the constant negative flux in is essentially a transient effect, 
which exists only in the period that the above stated condition holds. When the interaction time 77 exceeds 0(7"^), 
this approximation breaks down. To obtain the correct late-time behavior, one should take the limit rj^ ^ 00 before 
7 — > 0. Then ( AQ(t_)^ goes to zero. 

VIII. SUMMARY 

In this paper, we consider the Unruh-DeWitt detector theory in (3+1) dimensional spacetime. A uniformly accel- 
erated detector is modeled by a harmonic oscillator Q linearly coupled with a massless scalar field $. The cases with 
the coupling constant |Ao| less than the renormalized natural frequency \ flr\ of the detector are considered. 

We solved exactly the evolution equations for the combined system of a moving detector coupled to a quantum 
field in the Heisenberg picture, and from the evolution of the operators we can obtain complete information on the 
combined system. For the case that the initial state is a direct product of a coherent state for the detector and the 
Minkowski vacuum for the field, we worked out the exact two-point functions of the detector and similar functions 
of the field. By applying the coherent state for the detector, we can distinguish the classical behaviors from others. 
The quantum part of the coincidence limit of two-point functions, namely, the variances of Q and $, are determined 
by the detector and the field together. 

From the exact solutions, we were able to study the complete process from the initial transient to the final steady 
state. In particular, we can identify the time scales of transient behaviors analytically. When the coupling is turned 
on, the zero-point fluctuations of the free detector dissipates exponentially, then the vacuum fluctuations take over. 
The time scales for both processes are the same. Eventually the variance of Q saturates at a finite value, where 
the dissipation of the detector is balanced by the input from the vacuum fluctuations of the field. Even in the zero- 
acceleration limit, the variances of Q and Q, thus the ground state energy of the interacting detector, shifted from 
the ones for the free detector. This fluctuations- induced effect share a similar origin with that of the Lamb shift. 

The variance of Q yields an effective squared scalar charge, which induce a positive variance in the scalar field. This 
variance of the field contributes a positive radiated energy at the quantum level. However, the interference between 
the vacuum fluctuations and the retarded solution induced by the vacuum fluctuations screens part of the emission 
of quantum radiation. The time scale of the screening is proportional to 1/(7-1-0) for a > 7, where a is the proper 
acceleration and 7 is the damping constant proportional to A^. After the screening the renormalized Green's functions 
of the field are still non-zero in steady state. 

A quantum radiation formula determined at the null infinity of Minkowski space has been derived. We found that 
even in steady state there exists a positive radiated power of quantum nature emitted by the uniformly accelerated 
UD detector. For large a the radiated power is proportional to a'^Tjj, where Tjj is the Unruh temperature. This 
could be interpreted as a hint of the Unruh effect. However, the nearly constant negative flux obtained in Ref.,31j for 
(3+l)D case is essentially a transient effect. 

Only part of the radiation is connected to the internal energy of the detector. The total energy of the dressed 
detector and the radiated energy of a monopole radiation from the detector is conserved for every proper time interval 
after the coupling is turned on. The external agent which drives the detector's motion has no additional influence on 
this channel. Since the corresponding monopole radiation of quantum nature ceases in steady state, the hint of the 
Unruh effect in the late-time radiated power is not directly from the energy flux experienced by the detector in Unruh 
effect. This extends the result in Ref. |2Ct l27| that there is no emitted radiation of quantum origin in Unruh effect in 
(1-1-1) dimensional spacetime. 

Since all the relevant quantum and statistical information about the detector (atom) and the field can be obtained 
from the results presented here, when appropriately generalized, they are expected to be useful for addressing issues 
in atomic and optical schemes of quantum information processing, such as quantum decoherence, entanglement and 
teleportation. These investigations are in progress. 
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APPENDIX A: DERIVATIVES OF TWO-POINT FUNCTIONS OF THE FIELD 

From Gl^ in it is easy to see that 
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Note that the (5-functions at ry_ — 0, coming from the derivative of the step functions, have been neglected here. 

With + G^^ , one can write down in closed form of the interfering terms in the R- wedge of the Rindler space. 
Under the coincidence limit, it looks like, 
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with +e'"'+ and —6"''+ for a; in R and F-wedge, respectively. Note that is independent of 77+. Qxx is actually 
proportional to in H94() . Another observation is that, combining I|A1|I and (|A2|) . one finds that the divergent Ai 

term in Q is canceled by the one in ( Cp{rj_) )^. 

As for G**, the result is similar to (|Aip with (• • ■)^ being replaced by (• • •)^. Again it splits into the quantum and 
semiclassical parts as in Sec. IV B. 



APPENDIX B: A SIMPLER DERIVATION OF THE CONSERVATION LAW 



The conservation law (|109|) is directly obtained by arranging our somewhat complicated results. It can also be 
derived in a simpler way as follows. 

From the definition of the ground-state energy H82|) of the dressed detector, its first derivative of t is 



E{t) = moRe ( Q{t)Q{t) Q{t)Q{t) ) 



(Bl) 
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Introducing the equations of motion H38|l and H45|l to eliminate Q{t), one has 



E{t) ~ moRe <^ ( Q(r) 
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where the last term is a short-hand for 
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(There is an additional term AoQ(T)$in(y(T)) if a background field $in(a;) — ( <I'o(a;) )v is present.) Substituting ||1T)) 
and (|56(l into ljB3p . integrating out d'^fc with the help of Eq.ljHHl, then comparing with d^,dlG^°{x,x') in (ESJ where 
has the formal expression H93|l . one finds that (|B3|I is actually identical to —T%'^Q--{t)/ti. Hence 



E{T) = ~2^mo{Q\T)) 



(B4) 



Integrating both side from t = ti to t/, one ends up with Ea. (|109() . 
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